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We study the expansion of the surface thickness in the 2-dimensional lattice Sine 
Gordon model in powers of the fugacity z. Using the expansion to order z 2 , we derive 
lines of constant physics in the rough phase. We describe and test a VMR cluster 
algorithm for the Monte Carlo simulation of the model. The algorithm shows nearly 
no critical slowing down. We apply the algorithm in a comparison of our perturbative 
results with Monte Carlo data. 
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1 Introduction 



The 2- dimensional lattice Sine Gordon (SG) model is a prominent member of a large class 
of solid-on-solid (SOS) models that are believed to undergo a roughening transition of the 
Kosterlitz-Thouless type [1, 2]. The model can be considered as a Gaussian model with a 
perturbation that is periodic in the real variables ip x . The Hamiltonian is 

h ^ = Yb ^ (<p*-<Pv) 2 -v(<p). (i) 

1 <x,y> 

The sum runs over all nearest neighbour pairs < x,y > in the 2-dimensional lattice, and 

V(tp) = Z ^C0S(27T^). (2) 

X 

We have chosen units such that the Boltzmann factor is exp( — H). 

Compared to other SOS models 4 the SG model is very suitable for analytical calcula- 
tions, e.g., expansions in the fugacity z. 5 

The variables ip x have a natural interpretation as height variables. For sufficiently 
small temperature (3 a typical ip- configuration describes a more or less smooth surface. In 
the thermodynamic limit the surface gets localized and has a finite thickness. The model 
is in the smooth phase. 

For large /3, however, the surface can freely wander. Furthermore, the surface thickness 
squared diverges logarithmically when the area becomes large. The model is in the rough 
phase. 

The two phases are separated by a phase transition of infinite order, the roughening 
transition. The transition occurs at a critical line (3 c (z). 

To make things more precise, we define the surface thickness a. On a finite periodic 
square lattice A with N = L X L sites, let 

CT2 = ^2E<(^-^) 2 > = 2((^-<i>) 2 >. (3) 
$ is the average of ip over the entire lattice, 
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* = ^£^- (4) 

xEk 



A proof is presented in appendix 1 that for z = and L — > oo, 



<T 2 ^@-hiL + c, (5) 

7T 



4 For a Monte Carlo renormalization group investigation of several SOS models see ref. [3]. 
5 The name fugacity for the coupling constant z will become clear later when we consider the Coulomb 
gas representation of the Sine Gordon model, see also ref. [4]. 



1 



where c is some constant. 

The surface thickness can be used to distinguish between the two phases of the Sine 
Gordon model: In the limit L — > oo, for nonzero z, 

CT 2 J ^ilni + c', if /3 >&(*), (g) 
1 finite, else . 

This equation says that in the rough (large (3) phase, the long distance behaviour of the 
theory is that of a Gaussian model (z = 0) with infinite correlation length, but with a 
changed temperature /5 e ff that depends on the fugacity z. 

It is known that for z — > 0, the critical temperature goes to 2/ir, see, e.g., [4]. Accord- 
ingly, the critical line is given as the set of points (z,(3) with (3 e ff(z) = 2/ir. For every z, 
the critical point fi c {z) is the smallest value of (3 such that the surface thickness diverges 
for L — > oo. At (3 = (i c (z), the asymptotic ratio a 2 /ln(L) jumps from 2/ir 2 to zero. 
More generally, we shall define lines of constant physics in the rough phase of the model 
as those lines in the (z, /3)-plane that belong to the same value of /5 e ff . Models that lie on 
the same line have identical long distance behaviour (which is Gaussian in this particular 
model). 6 

In this paper, we shall derive an approximation of the lines of constant physics from 
the perturbative expansion of the surface thickness a to second order in z. 

It is always interesting to compare analytical results with numerical ones, especially 
when approximations are made in the analytical calculations. To make this comparison 
we employed a cluster algorithm for the simulation of the SG model. This algorithm is 
a new member in a class of cluster algorithms that we call VMR (valleys-to-mountains- 
reflection) algorithms [5]. In this paper, we shall describe the algorithm and study its 
dynamical critical properties. We use it to produce estimates for the surface thickness and 
other quantities that we can compare with the perturbative results. 

This paper is organized as follows. In section 2 we develop the perturbation theory 
for the surface thickness. Section 3 is devoted to the derivation of a formula for the lines 
of constant physics in the rough phase. We then turn in section 4 to the description 
and study of the VMR cluster algorithm for the SG model. In section 5 we compare the 
analytical results with those obtained with the Monte Carlo algorithm. 



2 Perturbation Theory for the Surface Thickness 

Expectation values (correlation functions) of observables in the Sine Gordon model are 
defined through 

/n(lo] , /rL^exp [-ff(y)] 0(y) 

{0(LP))= JlL^expt-i^)] • & 



6 According to D. Stauffer, one might call these lines also "isotrachs", a term composed from greek 
words, which means "lines of constant roughness". 
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The Hamiltonian H(<*p) is invariant under shifts ip x — f ip x + n for all x £ A, where n is an 
integer constant. As a consequence, expectation values are defined only for observables 
0{<~p) that have the same global symmetry: In this case the infinite contribution from the 
zero mode associated with the symmetry is exactly cancelled in eq. (7). 7 

For the perturbation theory to be done in this section, it is convenient to rewrite eq. (7) 
with the help of a Gaussian measure. An integration measure dupci'-p) lS defined through 
its "generating functional" 

d^cM exp[i(* )¥ 0] = { e Q XP t-W^*)] ; *Z* k ° = ° (8) 

The covariance C is a matrix with elements 

_ 1 e w(x-y) _ i 
xy ~ at 



N ^ p 2 



p = 4 — 2 cospi — 2 cosp2 i (9) 

where the pi, i = 1,2, are summed over the values {0,...,L — 1} • (2ir/L). The scalar 
product is defined by 

(<^) = ]>>^- ( 10 ) 

X 

It is shown in appendix 2 that the expectation value eq. (7) can be rewritten as 

[ W) fdfipciv) exp[T^)] • 1 ) 

The calculation of the interface thickness is started with the observation that for arbitrary 
observables 0((p), 



_ fdfipcjy) exp[F(y)] Q(<p) 
fdfj,p C ((p) exp[V(<p)] 

= fd$f d^cjf) S($ -jfJ2x <Px) exp[y(y)] Q{<p) 
J d<S> J d^ c (f) - j? Ex <Px) exp[y(¥»)] 

_ fd<S>f duprjy) exp[F($ + y)] + y) 
f d$ f dnf)r(<p) exp[F($ + v?)] 

Here, the Gaussian measure with covariance T is defined through 

d^M 0(<p) = lim / t7 ( 1'^ ) ° (y) ' ( 13 ) 

M^O J dfl VM ((f )S( J2x <Px) 



'12) 



r The Z-symmetry of the model is broken in the smooth phase, (6 < fl c (z). Here, in the thermodynamic 
limit, expectation values of observables that are not Z-invariant, are defined. 



3 



where vm = ( — A + M ) . It is not difficult to demonstrate that 

r I e w(x-y) 

T xy = J d[i r (<p) VxVy = E —~2 — • ( 14 ) 

Correlation functions as occuring in eq. (12) can be computed with the help of the gener- 
ating functional 

dfipriv) exp [i(k, (p)] = exp [-\j3{k, Tk)] . (15) 
The "effective observable" for the surface thickness is independent of 

+ = 2^. (16) 

We therefore have 

2 _ ? Idfipr(<p) <P 2 X fd$ exp[F($ + y)] 

fdiip r (<p) JdQexplV + ' 1 ) 

We expand the Boltzmann factor: 

z n 

exp[V($ + <p)] = ]T— E---E cos27r (^! +$)... cos 2^^ + $) (18) 



n>0 ^l,Sl X n ,S n 



(19) 



The "charges" s 4 - are summed over the values ±1. The <i>-integral leads to a ^-function 
that forces the sum of the s 4 - to be zero (neutrality condition). We shall indicate this 
constraint with a prime at the sums in the equations to follow. The Gaussian integrations 
in eq. (17) can now be done, leading to the representation 

V (^\ n 7' 

<x 2 = 2(3T 00 - K ^ w> ° * 4" • (20) 

Here, we have defined 

z 'n = ^ E • • • E ex p ( -/3 E - r A (E C °*A , (2i) 

■ x\,s\ x„,s„ y i<j y \i=i / 

and 

1 ' ' / ~ " \ 

z - = ^? E • • • E ex p -/^E^^ • (22) 



? For the definition of the fattice Lapfacian see appendix 2. 
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Furthermore, we have introduced the abbreviations 

/3 = (2vr) 2 /3 
23 2 

K = (2^' < 231 

In the derivation of eqs. (21) and (22) we exploited the fact that for neutral charge con- 
figurations s, 

n n 

Sj^XjXj s .j = 2 s iC XtXl Sj , (24) 
i,3 i<3 

and 

n n 

~] fpij s i = ^ ] C 0Xt Sj . (25) 

i i 

Eq. (20) is the Coulomb gas representation of the surface thickness. The gas is a neutral 
gas with integer charges ±1. A lattice site can be occupied by more than one charge. The 
representation is in the grand canonical ensemble, and z obviously plays the role of the 
fugacity for the charges S{. A small z means that the dominant contributions to the sum 
come from systems with only a few charges present (dilute Coulomb gas). Note that it 
follows from this representation that a finite fugacity always makes the surface thickness 
smaller, because only positive terms are subtracted from the z = result. 9 Furthermore, 
together with eq. (6) it follows that /5 e ff is smaller than (3. 

We now determine the coefficients in the expansion 

- 2 =Efir^- (26) 



One finds 



n>0 ^ 



a 2 = for n odd 

<?l = 2/3r 00 

a 2 = -KZ' 2 

a\ = —K(Z' A — Z' 2 Z2 



<J, 



6 — — (^6 — Z2Z4 — Z2Z4 + Z2Z2 ) 



a\ = ... (27) 
Let us explicitly write down Z' 2 , Z2, and Z' A : 



Z' 2 = J2 ex P (P C *v) ( C °* ~ C < 



\2 



J2 exp (fiC 0X ) J2(Co,,+y - C , y ) 2 , (28) 
x^o y 



3 We here assume that the infinite sum over n converges. This is probably the case in the rough phase. 
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Z 2 = L 2 J2 exp [PC 0X ) . (29) 

X 

Finally, 

^4 — ^ ^ ^ GXp | /3 (C x\X2 C X1X3 Cx^Xi C X2X3 C X 2Xi ~t~ Cx^Xi )| 

■ (Coa;i "I" C ox 2 ~ C 0Xs ~ C 0Xi ) (30) 



4 

^1,^2,^3,^4 



Of course, using translational invariance of C, the expression for Z' 4 can be transformed 
in a similar way as for n = 2. 

3 Lines of Constant Physics 

In the rough phase the behaviour of the 2- dimensional SG model at large distance is that 
of a Gaussian model (z = 0), but with an effective /5 e ff < fi. For the surface thickness this 
means that in the limit L — > oo, 



2 



lni + c'. (311 



We define lines of constant physics as the sets of points (z, (3) that belong to the same 
/3 e ff- Constant physics here means that the models lying on the same line have identical 
long distance behaviour. In this section we will use the expansion of the surface thickness 
to second order in z to determine the lines of constant physics. With eq. (28) and the 
observation that 

I e w x _ i 

J2( C o,x+y ~ C 0: y) 2 = -2— ^ /~ 2 ^ 2 (32) 



N (p 2 



one arrives at 



ft 2 i \ 1 p*P x — 1 

- 2 + -p (PC 0X ) - £ — ^ + 0(A . (33) 

^ > x^O p^O KF ' 

The first observation is that the 0{z 2 ) term of the z = surface thickness also behaves 
like lni for large L. To show this, we interchange the order of summation, 

/ \ 1 p ir P x — 1 1 1 / \ p ir P x — 1 

E -p (pcj) - E = ^ £ E -p (pcj) -— . (34) 

x^O p^O yI > p^O 1 x^O 1 
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Note that D p = A + 0(p 2 ), and A has an infinite volume limit: 

(3C 0X ) 

x^tO 1 

= -lJ2ex V ({3C 0X ) x 2 . (35) 

x^O 

The i-sum now runs over the infinite lattice. The 0(p 2 )-terms in D p will lead to corrections 
of a 2 that remain finite in the infinite volume limit. The limit p — > in eq. (35) is a bit 
subtle. Note that the limit of (px) 2 /p 2 alone strongly depends on the direction in that the 
limit is taken in the 2-dimensional p-space. However, the limit is to be performed under 
the sum over x. Now divide this sum into a sum over equivalence classes X of points x on 
that C ox takes the same value (these points are actually connected by lattice symmetries). 
Sum (px) 2 /p 2 over x £ X and then let p — > 0. Then the limit becomes independent of the 
direction in that the limit is taken, and one arrives at the result given in eq. (35). 
The large L approximation is thus 

° 2 a ° + z2 j£v h A °° + const + 0(z ' ) • (36) 



We apply the approximation formula [8 

1 

2~TT 



^--—(^1^1 + 1^2 + 7), (37) 



where 7 = 0.5772... denotes Euler's constant. The approximation for A^ then is 

Aoo « -iexp [-27r/3(|ln2 + 7)] ^l^l 2 " 2 ^ . (38) 

x^O 

s v ' 

E 2 (tv/3-1) 

E2 denotes the Epstein-^-function [9]. 

Note that replacing the "true" Coulomb propagator by its approximation eq. (37) will 
change A^ by a small finite amount (the approximation of the propagator has errors that 
decay like l/|a;|). However, A^ diverges for (3 — > 2 /it. Since the divergence comes from 
the large x contributions, the finite errors from the small distances can be neglected in 
this limit which will be of importance later. 
Recall that a 2 = 2(3T 00 . Therefore 

a 2 2T 00 (p - z 2 |(27r) 2 /3 2 exp [-2vr/3(§ In 2 + 7)] £ 2 (vr/3 - 1)) + const + 0(z 4 ) . (39) 

If we define a A/3 function through 

/3 eff = fi - z 2 ^fi(p) + 0(z 4 ) , (40) 
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we find 

A/3(/3) = i(27r) 2 /3 2 exp [-2vr/3(f In 2 + 7)] £ 2 (vr/3 - 1) . (41) 
o 

Let us now introduce a new variable 

2 

x = (3 - - . (42) 

7T 

x measures the distance from the critical point at vanishing fugacity z. For Ax(x) one 
obtains 



A; 



t(x) = i(2vr) 2 (1 + x^j exp (-2* (1 + x"j (§ In 2 + 7)) E 2 (l + ttx) . (43) 



We are now able to obtain the lines of constant physics. They will be parameterized in 
the form z(x, x eS ), where x eS = (3 eS - f: 

^ : '<«> = (t^)) 1 ' 2 ■ ,44 » 

In figure 1, we show the lines of constant physics for x e g = 0.0 . . .0.09. The case x e g = 0.0 
corresponds to the critical line. For the numerical evaluation of Ax(x) we used that [9] 
the Epstein (^-function E 2 , 

a( " = 5RT^' ,45) 

can be written as the product 

E 2 (s) = 4( R (s) P(s) , (46) 
where (r(s) denotes the Riemann ^-function, 



00 



( R (s)=J2n- s , (47) 



n=l 



and [10] 

00 1 poo / s_1 

«.)=E(-l)"(2»+l)-=jfw/ "^fy <«) 

Figure 1 looks like the famous flow diagram of Kosterlitz and Thouless. The strong 
similarity becomes more evident if we expand z(x, x e ff) around x = 0. Using the expansions 

Cr(s) = -^— + j + 0(s-1), (49) 
s — 1 

and [11] 

M = \(l + (7 + In 2vr + 2 In ( a - 1)) + O ((* - l) 2 ) , (50) 
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we find 

E 2 (l + irx) = ^ + k + 0(x), (51) 

with 

k = vr ^2 7 + In 2vr + 2 In = 2.5849817... . (52) 

If we plug in all our material, we find 

z(x,x eS ) = S ^Jx(x - x eS ) ^l + 0(x)j, (53) 

where S = exp((5/2) In 2 + 27) ~ 17.945. The critical line in the vicinity of the fixed point 
(z = 0, x = 0) is thus 

z(x,0) = Sx (l + 0(x)j . (54) 

For a comparison: The solutions of the Kosterlitz-Thouless equations in the rough phase 

are of the form 

Y = VE + X 2 , (55) 

with E < 0. Y and X are proportional to the fugacity and to (3 — 2/ir, respectively. 10 The 
case E = corresponds to the critical line. 

4 Cluster Algorithm for the Sine Gordon Model 

4.1 VMR for the Sine Gordon Model 

Cluster algorithms have proven successful in fighting critical slowing down (CSD) in the 
simulation of spin models. Cluster algorithms first occured in the pioneering work of 
Swendsen and Wang [6] in the context of the Ising model. 

Let us first describe the Swendsen- Wang algorithm as applied to an Ising model with 
partition function 

Z = ^2 ex P X! K <*,y> a x a y ' ( 56 ) 
<r x =±i \<x,y> j 

where <x,y> is the bond connecting the sites x and y. 

In the first step of the algorithm all bonds of the lattice are deleted or frozen. A bond 
< x,y > is deleted with probability 



1 exp(-K <Xiy> (l + <T x <T y )) if K <XiV> >0, ( 
Pto> " 1 exp(+A<^ > (1 - o x o y )) if K <Xiy> <0, [b<) 



10 Note, however, that the corresponding constants are nonuniversal, i.e. depending on the cutoff scheme. 
They can not be determined by a triviaf computation. 
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or frozen with probability Pf< x ,y> = 1 — Pd<x,y>- In a second step, the clusters of sites 
connected by frozen bonds are identified. Finally, each cluster is flipped with probability 
one half. 

Wolff proposed the so called single cluster version of this algorithm [7]. Here in one 
update step only one cluster is built. One randomly selects a seed x for this cluster 
and connects neighbouring sites with the above given probability pf to this site. This 
procedure is continued with the neighbours of the new sites until no further sites join the 
cluster. Then one flips the spins in this cluster with probability one. 

Cluster algorithms for the simulations of SOS models models were proposed in [5]. The 
idea for these algorithms came from the picture of SOS configurations as landscapes with 
hills and valleys. Large scale changes of a configuration are done by choosing a horizontal 
reflection plane, considering the connected regions above the plane (hills) and below the 
plane (valleys), and reflecting them through the plane independently, with an appropriate 
probability. 

We shall describe the details of the algorithm in the language of the Discrete Gaussian 
(DG) model. The modifications necessary for the Sine Gordon model will be given later. 
Note that the DG model might be considered as the z — > oo limit of the SG model. The 
DG model has the partition function 



where h x is an integer, and the sum runs over nearest neighbour pairs. 

Let us denote the height of the horizontal reflection plane by M. A reflection of h. 
with respect to M means 



Obviously, M has to be either an integer or a half-integer. The reflection can be expressed 
in terms of embedded Ising variables, s x = ±1, defined by the decomposition 



We get the partition function of the embedded Ising model by inserting eq. (60) in eq. (58), 




(58) 




(59) 



h x = s x \h x - M\ + M . 



(60) 




(61) 



with 



K. 



<x,y> 



= 2K 



h x - M\ \h y - M\ . 



(62) 



Hence we get from eq. (57) 



Pd<x,y> 



exp (-2K DG \h x - M\ \h y - M\ (1 - s x s y ) ) j . 



(63) 
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It is important that the reflection plane lies with reasonable probability within the vertical 
bounds of the SOS surface. In [5], we introduced two different types of algorithms to ensure 
this. 11 In the H-algorithm one chooses M in the neighbourhood of the height of the seed 
spin h Xo . More precisely, M is selected with a priori probability 

P < M > = { I ^ M = h " ±h < 64 > 

The H-algorithm alone is not successful in eliminating CSD. The situation is improved by 
including reflections with respect to an integer reflection plane M. In the I-algorithm, one 
selects an integer reflection plane M by simply setting it equal to the height variable at a 
randomly chosen site not equal to the seed: 

M = h yo , y ol L x . (65) 

Since the I-algorithm is not ergodic, we combined it with the H-algorithm. For the com- 
bined algorithm we did not find any CSD in the simulation of the DG model both in the 
rough and in the smooth phase. 

We now turn to the SG model. The Hamiltonian is 

= ;A (<Px ~ Vv? ~ zY, cos ( 2irL P^) ■ ( 66 ) 

' <^,y> x 
We use the embedding of Ising variables 

<Px= s x \<p x -M\ + M. (67) 

Because of its symmetry properties the SG potential remains unchanged if we perform 
cluster flips with respect to a reflection plane M that is either integer or half integer. Our 
version of the H-algorithm for the SG model is specified by 

™={ eke M = ™' tW " ]± ^ (68) 

and for the I-algorithm by 

M = nmt((p yo ) , y ^ x , (69) 

where y is again a randomly chosen site not equal to the seed x . Furthermore, nint^^) 
is the integer n such that \n — <~p x \ is minimal. The delete probabilities are now given by 

Pd<x, y > = exp ^-i \h x - M\ \h y - M\ (1 - s x s y )^j . (70) 

The algorithm described so far is not ergodic. We therefore include in our update procedure 
a local Metropolis updating of the spins. 



For a discussion of the detailed balance condition see the same reference. 
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4.2 Performance of the Algorithm 

We made an attempt to estimate the dynamical critical exponent of our VMR algorithm 
in the rough phase. To this end we measured various quantities. Local observables were 

E = ((ip x — (Py) 2 ) , with x,y nearest neighbours (71) 

and 

C = (cos(27r^)> . (72) 

As nonlocal quantities we measured the surface thickness a and a set of "block spin 
observables". To this end, the lattice was divided into four blocks x' of size L/2 X L/2. 
Block spins <f> x i were then defined as the averages of p over these blocks. We looked at 12 

A ^ = \ £ ((^'-<V) 2 )> (73) 

<x',y'> 

and 

^3,2 = J ]T (C0S(27T<M) • (74) 
x< 

The lattice average $ of <p was already defined in the introduction. We also measured the 
quantity 

A 3 ,i = (cos(2tt$)) . (75) 

The efficiency of a stochastic algorithm can be characterized by the integrated autocorre- 
lation time 

-i oo 

r = g £ P(<) > (76) 

t= — CO 

where the normalized autocorrelation function p(t) of an observable is given by 

Pit) - (Q.-Q^)-(Q) 2 (77) 

p( ) ~ (0 2 )-(0) 2 ' ( ' 

We calculated the integrated autocorrelation times r with a self-consistent truncation 
window of width 6r max for the observables specified above, where r max was the largest 
autocorrelation time found in the set of quantities considered. 

In the following we shall discuss the performance of the algorithm at (3 = 0.75 and 
z = up to z = 1.5 in steps of 0.5. We made runs on lattices from L = 8 to L = 256. We 
always made 100000 measurements of the quantities specified above, the measurements 
separated by an update step that consisted of the generation and flip of one H-cluster, one 
I-cluster, and one Metropolis sweep. The proposed change of the spins in the Metropolis 
step was always taken with uniform probability from the interval ( — 0.5,0.5). 



2 The notations are chosen to be consistent with the notations of ref. [3]. 
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It turned out that the surface thickness and A^^ had the largest autocorrelations. 
The integrated autocorrelation times for these quantities are given in table 1. Another 
observation that might be of general interest: The autocorrelation function of A^^ and 
the autocorrelation functions of the other block observables showed a clean exponential 
decay already for moderately small /. For the other quantities we could not observe such 
a clean exponential decay. 

For the ^-values considered, the point (3 = 0.75 is definitely in the rough phase, i.e. 
the correlation length is infinite. The r's are therefore expected to scale with the lattice 
size L like r oc L z , where z is the dynamical critical exponent (not to be confused with 
the fugacity). 

We fitted our data to this law. When the smallest lattice sizes were discarded, we 
obtained fits with good x 2 /dof. Our results are given in table 2. The results for the 
different fugacities are fairly consistent with each other. Furthermore, an exponent of 
order 0.3 indicates that the CSD has been drastically reduced. 

Let us now give a few more details on the runs at (3 = 0.75 and z = 1.0. The clusters 
built with the H-algorithm had an average size (in units of L 2 ) of 0.51 for L = 8 to 0.33 
for L = 256. For the clusters of the I-algorithm the average size was fairly independent 
of the lattice size, always close to 0.32. The acceptance rate in the Metropolis step was 
approximately 60 per cent throughout. 

A word on the resources needed for this study: the runs with L = 256 needed something 
like 23 hours on an HP9000/735 workstation. 

We did not attempt to estimate the dynamical critical exponent in the smooth phase, 
but we have no doubt that the performance is comparably good there. 

5 Comparison of Perturbation Theory and Monte Carlo 
Results 

5.1 Accuracy of the Perturbation Theory 

We checked the accuracy of eq. (33) by a direct comparison with Monte Carlo results for 
the surface thickness obtained with the VMR cluster algorithm described above. 

In table 1 we compare the results at (3 = 0.75. For z < 1 the deviations are fairly 
small. For z = 1.5 the relative error of the perturbative result is of order 1.5 per cent. 

A similar statement is true for a set of (3 values that we chose to be close to the 
roughening point (3 c (z) (compare subsection 5.2). Table 3 shows our results for the points 
(z,f3) = (0.5,0.665), (1.0,0.700), and (1.5,0.720). 

Table 4 shows the same comparison for three (3 values in the smooth phase. As is to be 
expected from the fact that <7q and <rf diverge logarithmically with L, the approximation 
fails completely for large L: the breaking of the symmetry under global integer shifts and 
the finiteness of a 2 for L — > oo are not reproduced by the perturbation theory. 
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5.2 The Critical Line 



We tried to determine fi c {z) for z = 0.5, 1.0 and 1.5 using the matching method that we 
introduced in [3]. This method is based on a comparison of the blocked observables of the 
model in question with that of the BCSOS model. The BCSOS model is exactly solved and 
known to undergo a KT phase transition. The basic idea of our matching procedure is to 
simulate the BCSOS model at its known roughening temperature. One measures various 
block observables, e.g., the Aij introduced in subsection 4.2. One then searches for the 
critical coupling (3 C and a scale transformation of the SOS model under consideration such 
that the blocked SOS observables match with that of the critical BCSOS model. 

It was demonstrated in ref. [3], that the matching with the BCSOS model can be per- 
formed by using just two block observables, e.g., an( i ^-3,2- The matching conditions 
then fix (3 C and the scale transformation. The matching of all the other block observables 
is then a consequence of universality and can be demonstrated. 

We tried to use the block BCSOS observables reported in [3] for a matching analysis 
of the SG model. We restricted ourselves to the quantities an( i ^-3,2- 

It turned out that the values of A^^ f° r the SG model close to the expected critical 
temperature are smaller than those of the BCSOS model on lattices of size up to L = 128 
[3]. Since A^^ vanishes like the inverse of the logarithm of the lattice size there is no hope 
to simulate the BCSOS model on lattices so large that the direct matching with the SG 
model at the z values under consideration can be performed. 

From KT theory one expects that to leading order A\^ is a linear function of A^^- In 
figure 2 we plot for the BCSOS model the quantity D\p as a function of A^^, where 

DiM = t' 2 mh ^ 2 W • (78) 

In [3] we demonstrated that replacing A\^ by -Di,2 reduces finite lattice-spacing artefacts 
considerably. In addition to the BCSOS results the z = result is known. It is given by 
As t 2 = and A\^ = 0.075425. It seems that the linear interpolation between this z = 
point and the L = 128 point of the BCSOS model is a reasonable approximation to the 
curve. On obtains A\^ = 0.075425 + 0.034(2)^3^. fi c is now computed from the condition 
that Ai^ift) and A^^ifi) are elements of the curve. 

We simulated the SG model at /3's close to the expected value for (3 C on lattices up 
to L = 64. We computed the values for the two observables in the neighbourhood of 
this (3 using the reweighting technique. Instead of A\^ itself we again considered D\p m 
order to reduce finite lattice-spacing artefacts. The results for (3 C stemming from different 
lattice sizes are summarized in table 5. We give two error bars. The first stems from the 
statistical error of the SG data itself, while the second is due to the error bar of the slope 
of the critical ^^2(^3^) curve. Even the result from L = 8 is consistent with the results 
from larger lattice sizes within the statistical errors. 

From the exact result at z = and the estimate /3 C (0.5) = 0.670(2) we derive a slope 
of the critical line of S = 14.98(0.9). This is obviously inconsistent with our estimate 
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5 = 17.945 from the second order perturbation theory. 

Beyond a three standard deviation statistical error there are two possible sources for 
a systematic error. First, the linear approximation of ^1^(^-3, 2) might be inadequate. 
Second, we were not able to check whether higher orders in z of the surface thickness 
expansion can modify the slope of the critical curve. Such a modification would happen 
if the analogue of at order z n would diverge like l/a;( n_1 ) for x — > 0. 

6 Summary and Conclusion 

We have presented a computation of the lines of constant physics of the SG model in 
the rough phase. The calculation was based on an expansion of the surface thickness to 
second order in the fugacity z. The comparison of the surface thickness data with the 
ones obtained with a very efficient VMR cluster algorithm revealed the reliability of the 
expansion in the small fugacity regime. When comparing the slope of the critical line with 
the result obtained from a RG matching with the BCSOS model, we found a significant 
deviation. We cannot presently decide whether this effect is due to the higher order terms 
in z, or whether there is another, yet undiscovered reason for the deviation. 
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Appendix 1: The z = Surface Thickness 

The surface thickness for z = is given by 

* 2 = ^£<(^-<AOV- (79) 

X 

It is not difficult to demonstrate that 

°o = wjf E 4 , (80) 

where p 2 = 4 — 2cospi — 2cosp2 5 an( i the Pi, i = 1,2, are summed over the values 
{0,...,Z-1}-(2tt/Z). 

We want to prove the following 

Statement: For L — > oo, the surface thickness a 2 goes like 

3 

a* -> - hiL + c. (81) 
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For the proof of the statement we need the following definitions (for formal reasons let us 
assume that L is even): 



Al = jP = (Pi,P2) - Pi = n t — , -(- - 1) < n t < - 1 



(82) 



(83) 



and 

A~ = jp = (pi,P2) : Pi = ra *"~j~> _00 ^ ra * ^ °°} • 
Now we define the quantities 

Q(L) = jl E h (84) 
P eA i \{o} i ' 

and 

2 2 

R(L,a) = ± J2 (85) 

p£A-\{0} ^ 

Note that a provides an ultraviolet cut-off in the second sum. To prove the statement we 
need the following 



Lemma: For a > 0, we have 



lim {Q{L) - R(L, a)} < oo . (86) 

L— >oo 



Proof of the Lemma: We write the difference of Q and R as a sum of two contributions, 



1 ^ f 1 11 1 ^ e" aV 



Q(L)-R(L,a) = - £ --- -- ^ _. (87) 
^ P eA L \{o} lp p J ^ pga^\Aj, ^ 

' v ' v 

In the limit L — > oo, the first term becomes 

r d 2 p f 1 1 1 

because 

lim | — - —1 = finite. (89) 
p^o [p 2 p 2 ) 

Let us define 

G := {p : \ Pi \ > vr} . (90) 
In the limit of large L the second term becomes 

Jg (27r) 2 p 2 27r A q 
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It follows from the Lemma that the parts of Q(L) and R(L, a) that diverge when L — > oo 
are identical (since their difference stays finite in this limit). Employing the standard 
"Feynman" trick, we get 

R(L,a) = — }2 —T- 

p£A-\{0} P 

= Tlf^E e >V -l}^ ( 92 ) 



Here, A denotes the Jacobi-^-function, 



A(s):=J2 e ~ vn S ■ (93) 

neZ 



After two simple substitutions we find 



1 f°° 

^' a)= 4^yw^ w)2 - i} (94) 

i 2 



For / — > we have 



A(t)^^{l + e-T + ...}. (95) 
It is therefore instructive to split the integral representation of R as follows: 

R(L,a) = -Ljl^dtiAitf-]} 

L 2 

1 fOO 

+ -/ 

If 1 1 

+ ~T 2 dt il ~ !>• (96) 

L 2 

The first two terms stay finite when L — > oo, the third integral diverges logarithmically: 

a) finite + — InL. (97) 

27T 

This is what we wanted to show. 



Appendix 2: Gaussian Measures 

Expectation values in the Sine Gordon model are defined by 



{oi } = lim /n^expt-g^y) _ (98) 
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Here, we have introduced a symmetry breaking "mass" term that serves temporarily as a 
regulator of the zero mode: 



1 M 2 
Hm(^) = ^ J2 (<P* ~ ^yf + — Z>* " ■ 

' < X ,y> x 



We define the scalar product 
and a lattice Laplacian A, 



(A^ = ifv ~ fx) , 

y.nn.x 



(99) 



:ioo) 



:ion 



where the sum runs over the sites y that are nearest neighbours of x. Eq. (99) can now 
be rewritten as 



A 

7 



f) ~ V(<p) • 



:io2) 



Now we write 



(0(<p)) = lim 

JTL^exp [-i(^,^V) + ^(^)J / JTL^exp [-i(^,^V)J 



Here, we have introduced the abbreviation 



A 



:io3) 



(104) 



The matrix j?m has an inverse vm that can be explicitly calculated by Fourier transfor- 
mation: 

1 e w(x-y) 

Vm ^=n\ p/p + M*- (105) 

Here, the pi, i = 1,2, run over the values {0, 1} • (2ir/L), and 

p 2 = 4 — 2 cospi — 2 cosp2 • 
The Gaussian measure with covariance vm is defined through 

JTL^exp \-\{ip,v^Lp) \ 0(<p) 



:io6) 



With this definition we can write 



/ iL^exp l-UfjVjJtp) 



(107) 



(C%)) = lim 



fd» VM (<p) exp c%) 



M-,o fdfj, VM ((p)ex.p[V((p)] 



:io8) 
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The generating functional of this measure can be explicitly calculated: 

dfj, VM ((p) exp [i(k, <p)] = exp [~\{k, v M k)] . 



:io9) 



All correlation functions can be derived from this relation. Let us now see what happens 
in the limit M 0. We split 



vm. 



xy 



1 e ip(x-v) _1 1 



1 



N 



p 2 /f3 + M 2 N^p 2 /f3 + M 2 



:no) 



The first term stays finite when M — > 0, the second (to be called K(M)) diverges to +oo. 
Therefore 



dnpc{f) exp [i(k, <p)\ = Jim o / dfj, VM (<p) exp [i(k, <p)] 



exp [—^(3(k,Ck)] lim exp 



■\k(m) {y^k x 

\ x 



exp[-^(k,Ck)] ,if£ x k x = 
, else . 



'1111 



The matrix elements of C are given by 

C 



'xy 



1 e ip{x-y) _ I 
~N ^ 



p2 



References 

[1] To give only a few references: 

J.M. Kosterlitz and D.J. Thouless, J. Phys. C6, 1181 (1973); 

J.M. Kosterlitz, J. Phys. C7, 1046 (1974); 

S.T. Chui and J.D. Weeks, Phys. Rev. B14, 4978 (1976); 

J.V. Jose, L.P. Kadanoff, S. Kirkpatrick, and D.R. Nelson, 
Phys. Rev. B16, 1217 (1977); 

T. Ohta and K. Kawasaki, Prog. Theor. Phys. 60, 365 (1978); 

D.J. Amit, Y.Y. Goldschmidt, and G. Grinstein, J. Phys. A13, 585 (1980). 

[2] D.B. Abraham, 'Surface Structures and Phase Transitions - Exact Results', in: Phase 
Transitions and Critical Phenomena, Vol. 10, C. Domb and J.L. Lebowitz, eds., Aca- 
demic 1986; 

H. van Beijeren and I. Nolden, 'The Roughening Transition', in: Topics in Current 
Physics, Vol. 43, W. Schommers and P. van Blankenhagen, eds., Springer 1987. 



19 



[3] M. Hasenbusch, M. Marcu, and K. Pinn, Physica A 208, 124 (1994). 
[4] S. Samuel, Phys. Rev. D 18, 1916 (1978). 

[5] H.G. Evertz, M. Hasenbusch, M. Marcu, K. Pinn, and S. Solomon, 
Phys. Lett. B254, 185 (1991); 

H.G. Evertz, M. Hasenbusch, G. Lana, M. Marcu, and K. Pinn, 
Phys. Rev. B 46, 10472 (1992). 

[6] R.H. Swendsen and J.S. Wang, Phys. Rev. Lett. 58, 86 (1987). 

[7] U. Wolff, Phys. Rev. Lett. 62, 361 (1989). 

[8] C. Itzykson, J.M. Drouffe, 'Statistical Field Theory', Cambridge University Press 
1989. 

[9] G.H. Hardy, Mess. Math. 49, 85 (1919). 

[10] W. Magnus, F. Oberhettinger, F.G. Tricomi, 'Higher Transcendental Functions', 
Vol. 1, A. Erdelyi, ed., McGraw-Hill, New York, 1953. 

[11] G. Miinster, private communication. 



20 



z 


7" 

L 


-43,1 






2 


2 

CTp T 


0.0 


8 


0.0087(57) 


3.27(07) 


1.37(4) 


0.56751(75) 


0.568942 


0.0 


16 


-0.0085(70) 


4.78(12) 


1.75(5) 


0.73483(81) 


0.734887 


0.0 


32 


0.0000(81) 


6.56(18) 


2.41(11) 


0.89919(92) 


0.900489 


0.0 


64 


-0.0005(87) 


7.53(24) 


2.72(13) 


1.06513(98) 


1.065998 


0.0 


128 


0.007(10) 


9.71(41) 


3.51(18) 


1.2324(11) 


-i O O "1 /I o o 

1.231483 


0.0 


256 


0.009(11) 


-1-1 r" n / *~» /-» \ 

11.58(39) 


3.78(21) 


"1 / "1 

1.3972(12) 


1.396961 


0.5 


8 


0.0571(58) 


3.38(08) 


1.45(5) 


0.56602(78) 


0.566367 


0.5 


16 


0.0539(70) 


A C\ O / -1 O \ 

4.83(13) 


1.73(5) 


0.73104(80) 


0.731812 


0.5 


32 


0.0258(79) 


6.19(19) 


2.31(10) 


0.89795(92) 


0.896745 


0.5 


64 


0.0307(86) 


7.35(22) 


2.63(11) 


1.05992(95) 


-i r\ f -i /i i~) /i 

1.061484 


0.5 


128 


0.0162(97) 


9.25(38) 


3.10(16) 


1.2264(10) 


1.226138 


0.5 


256 


n aat / 1 a\ 

0.005(10) 


10.67(41) 


3.89(17) 


-i o AO r / 1 o\ 

1.3925(12) 


-1 o AA "7 1 O 

1.390748 


1.0 


8 


0.1021(57) 


3.30(07) 


1.61(4) 


0.56079(85) 


0.558644 


1.0 


16 


0.0849(69) 


4.69(10) 


1.96(6) 


0.72417(87) 


0.722587 


1.0 


32 


0.0694(78) 


6.05(22) 


2.31(7) 


0.88976(92) 


0.896745 


1.0 


64 


0.0518(89) 


7.91(28) 


2.69(12) 


1.05077(96) 


1.061484 


1.0 


128 


0.0370(91) 


8.31(33) 


3.42(16) 


1.2148(11) 


1.226138 


1.0 


256 


0.023(11) 


11.29(39) 


3.98(19) 


1.3776(12) 


1.372110 


1.5 


8 


0.1466(57) 


3.34(08) 


1.62(4) 


0.55412(86) 


0.545772 


1.5 


16 


0.1114(69) 


4.79(13) 


1.92(7) 


0.71853(86) 


0.707212 


1.5 


32 


0.0984(77) 


5.91(17) 


2.36(8) 


0.87879(92) 


0.866793 


1.5 


64 


0.0779(90) 


8.02(24) 


2.73(13) 


1.04003(97) 


1.025370 


1.5 


128 


0.0786(97) 


9.28(27) 


3.69(20) 


1.2015(11) 


1.183376 


1.5 


256 


0.0634(99) 


10.00(41) 


3.82(20) 


1.3594(11) 


1.341046 



Table 1: Monte Carlo estimates for the quantities and a 2 for different 
lattice sizes L and fugacities z. (3 is 0.75 always. The Monte Carlo results 
for the surface thickness squared can be compared with leading order per- 
turbation theory (PT) given in the last column. The r's give the integrated 
autocorrelation times. 



z 


z A 3 ,i 


Z a 2 


0.0 


0.282(18) 


0.237(28) 


0.5 


0.290(15) 


0.285(17) 


1.0 


0.310(13) 


0.257(18) 


1.5 


0.294(15) 


0.264(20) 



Table 2: Estimates for the dynamical critical exponents as obtained from the 
integrated autocorrelation times of the quantities and a 2 . 
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Table 3: Monte Carlo estimates for the surface thickness close to (i c (z), and 
comparison with perturbation theory. 
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Table 4: Monte Carlo estimates for the surface thickness in the smooth phase, 
and comparison with perturbation theory. 
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Table 5: Estimates for the critical couplings fi c {z) from matching of an( i 
^3,2- Pc f° r ^ = oo is the roughening coupling for the Discrete Gaussian 
model [3]. 
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Figure 1: Lines of constant physics in the rough phase of the 2- dimensional lattice Sine 
Gordon model for x e g = 0.0 . . .0.09. 
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BCSOS at K c 




Figure 2: The quantity D\p as a function of A^^ f° r the BCSOS model at criticality. The 
data are taken from ref. [3]. 
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